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Abstract 



We critically examine the validity of the Weizsacker- Williams approxima- 
tion in electron-hadron collisions. We show that in its commonly used form it 
can lead to large errors, and we show how to improve it in order to get accurate 
results. In particular, we present an improved form that is valid beyond the 
leading logarithmic approximation in the case when a small-angle cut is applied 
to the scattered electron. Furthermore we include comparisons of the approx- 
imate expressions with the exact electroproduction calculation in the case of 
heavy-quark production. 
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Electroproduction phenomena are usually computed by assuming that the incom- 
ing electron beam can be considered to be equivalent to a photon broad-band beam 
PI. The photon distribution in the electron is given by the formula 



[■(e)/ - 



27r 



y mjy^ 



where 5* is the electron-hadron centre-of-mass energy squared, Smin is the minimum 
value for the invariant mass squared of the produced hadronic system, me is the 
electron mass, and y is the fraction of the electron energy carried by the photon. 
Several examples of applications of this formula, as well as discussions of its range 
of validity, can be found in ref. 0. We have explicitly indicated in eq. (|l]) that the 
neglected terms are of order 1 with respect to the logarithmically-enhanced one. The 
logarithmic term is typically of order 15 to 20 at HERA energies, and therefore one 
expects the error associated with this approximation to be of the order of 5 to 10%. 
It was noted in ref. that the error can in fact be much larger than this, unless one 
takes appropriately into account the dynamical effects that modify the "effective" 
upper scale entering the logarithm in eq. (^. 

In this paper, we show that several modifications to eq. (|1|) are needed in order 
to meet the needs of the photoproduction experimental conditions one encounters at 
HERA. There a photoproduction event is usually defined with an appropriate anti-tag 
condition (i.e. outgoing electrons above a given angular cut are vetoed). Therefore, we 
reconsider the derivation of the Weizsacker- Williams distribution, taking into account 
the particular experimental conditions of HERA. First of all, we derive a Weizsacker- 
Williams formula for the case when a small-angle cut is applied to the outgoing 
electron, which is valid beyond the leading logarithmic approximation. 

Let us consider the electroproduction process 

e{p)+p{k)^e{p')+X, (2) 

where p is a massless parton (fc^ = 0), X is a generic hadronic system, and = p'^ = 
rrtg. The cross section for the process is given by 

8k -p g4 (27r)32E' ' ^ ^ 

where 

q = p — p'. (4) 
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The tensor T^j, is the electron tensor, given by 

T^u = 4 {^q^g^,u + P^,p'y + PuP'^^j ■ (5) 

Notice that 

= g^T^, = 0, (6) 

since p ■ q = —p' ■ q = q'^/2. 

The tensor W'^'^ is the hadron tensor. Exploiting current conservation, which 
requires q^jW^^ = q^W^^ = 0, it can be decomposed as 

= W,iq\ k.q) + ^) - k.q) - ^g'^) - ^g^) 

(7) 

(a term proportional to e^yp^q'^k'^ has not been included, since it gives zero contribu- 
tion in this case). In the limit g^ — 0, W^^ must be an analytic function of g^. By 
requiring that q^W^^ vanish for g^ = 0, we obtain 

W2{q\k-q) = Wi{Q,k-q) + 0{q^). (8) 

We drop for the moment the 0(g^) terms. We will discuss later their effect, which 
in general leads to non-universal and non-logarithmic corrections to the factorised 
Weizsacker- Williams approximation. From eqs. (|]), and (|^) we find 



1^'^%, = -4W^i(0,A;-g) 

where 



2ml + q^ 



y2 



(9) 



k ■ q k ■ p' , ^ 

Notice that y is precisely the fraction z of longitudinal momentum carried by the 
photon. In fact, ii k = 0, 0, — 1) and q = zp + q, with k ■ q = one finds 

y={k- q)/{k-p) = z. 

We must now find a convenient expression for the phase space of the scattered 
electron. We have 

^ = 271 (3'^ E' dp' d cose, (11) 
where the trivial azimuth integration has already been carried out, and 



p' = {E', 0, E'I3' sin 9, E'(3' cos 6) 



(12) 
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p={E, 0,0, E(3) 



1-^ B 



E^' 



(13) 
(14) 



We now trade the two variables {p', cos 6) for (g^, y). From the definitions given above, 
we find 



g2 = 2ml - (1-/3/3' cos ( 

y 



E'(l + /?'cos( 



(15) 
(16) 



and it is easy to prove that the Jacobian of this change of variables is simply 2E /3 
Therefore 



= ndq dy. 

We are now ready to compute the cross section in eq. (0). We get 

aemWi{0,k-q) 



dap 



^ l + {l-yf 



27r Ak-p I 
where, as usual, = e^/(47r). Integration in gives 

dcTep = a^p{q, k)f^''\y)dy, 

where 



dq dy. 



27r 



2mly 



et> \ 2 2 



and 



a^p{q,k) 



1 + (1 -y)2 o2 

y ^max 

g^uW^" _ WiiO,k-q) 



8k ■ q Ak ■ q 

is the cross section for the process 7(g) + p{k) X for an on-shell photon. 



(17) 

(18) 

(19) 
(20) 
(21) 



We must determine the two integration bounds q^ax g^j„. For 6* ^ 1 we have. 



from eq. 



where 



E' 



+ 



2\2 

eJ n2 



2A ' 8A^ 
A = E{l + f3){l-y) 



and eq. (0) becomes 



mly"^ E{1 + (3){A^ - m\ 
l-y 4^3 



2\2 
e) q2 



(22) 
(23) 

(24) 
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The value of g^„^ is obtained by taking 9 = 0, namely 



^max 



2 2 



(25) 



Analogously, the value of is obtained when 9 equals its maximum value 9c- If 



6'c ^ 1 we can use eq. ( p4D to obtain 

E{1 + P){A^ - mlY 



2 2 
2 2 

m^y 



4^3 



+ 



E\l-y)9l + 0{E'9t,ml9l,mt/E 



(26) 



The function f!f\y) takes the form 



Or 



27r 
1 + (1 



+ 



2(1 -y) 



mly 



E'^{l-yf9l + mly'^ 
E\l - y f9l + mly'' 



y 



log 



mly'' 



+ 0{9lml/E' 



(27) 



Observe that the uncertainty is small with respect to 1. Also notice that the non- 
logarithmic term is singular in y and therefore represents a non-negligible correction. 
This subleading term is universal, in the sense that it is independent of the particular 
hard process we are considering. 

We now comment on the effect of the 0{q'') terms in eq. (§). It is straightforward 
to verify that these terms will add a contribution to the differential cross section given 
by 

.n=o \q-k) q-k 

where the coefficients A'^"'\q ■ k) are obtained from the derivatives of the W2 form 
factor evaluated at = 0, and are of order 1. For simplicity we neglected terms that 
are finite for y ^ 0, proportional to derivatives of both Wi and W2-, and for which 
the following conclusions hold as well. Upon integration over g^, we are left with 



Ma, 



1 



ep 



27r 2k -q 



dy 

y 



(28) 



Ma, 



ep 



27r 2k -q 



E 

Ln=l 



-1) 



n 



dy 

y 



(29) 



where we neglected the contributions proportional to qlnax^ which are suppressed by 
powers of ml/s. In the case of a small angular cut we have 



\q'^ ■ I 



< 



E'9l 



2q ■ k S, 



(30) 



-5- 



For example, in the case of charm production and using E = 30 GeV, 6c = 5x 10^^ we 
would have l^minl ■ k) < 2.5 x 10~^. Therefore these non-factorisable corrections 
are negligible, at least formally, with respect to the factorisable ones. 

If no angular cut is applied, however, the ratio \q'^in \ /(2g ■ k), although limited 
from above by 1, can be of order 1. In this case, non-factorisable corrections to the 
Weizsacker- Williams approximation are of order 1/y. 

As an application of the formulae discussed above, we consider the case of heavy- 
quark electroproduction. The total cross section for this process has recently been 
calculated at next-to- leading order in QCD by E. Laenen et al. in ref. [Q. For 
simplicity we will perform our comparisons between exact and approximated results 
at the leading order in as- For our goals this is not a limitation, as higher-order QCD 
corrections should not affect the properties of the photon density inside the electron 
and therefore will not change our results qualitatively. 

In table |I] we show the total production cross section of charm and bottom quark 
pairs at HERA {E = 30 GeV, ^/S = 314 GeV), with the apphcation of different 
angular cuts. We used rUc = 1.5 GeV, = 4.75 GeV, and the MRSDO parton 
distributions ^ for the proton. The matrix elements for the exact leading-order 
calculation are taken from ref. IQ. 

As can be seen from table |l|, the presence of the non-logarithmic term improves 
the agreement between the exact result and the Weizsacker- Williams approximation 
from the level of 6-7% to the level of 1% in the case of charm, and from 5-7% to 
0.1% the case of the bottom (apart from the case 6c = 0.5, which is too large for the 
small-angle approximation to work properly). 

Table |^ is similar, but evaluated for a configuration with electron beam energy 
E = 100 GeV and v^=l TeV. We see that with the same angular cuts but higher 
energy, the approximation becomes worse (this can easily be explained from eq. (pOD). 

Let us now consider the case when no angular cut is imposed on the scattered 
electron, and the angular integration is performed over the whole phase space. In 
this case, the minimum value of the photon virtuality can be computed in a 
simple way by imposing that the invariant mass of the produced hadronic system be 
bounded from below: 

(g + kf > Smin, (31) 

where, for example, Smin = 4M^ in the case of heavy-quark pair production. Equa- 
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Charm 




(7exact (nb) 


a^w only log (nb) 


(%) 


aww (nb) 


(%) 


5 X 10-1 


353.3 


420.7 


19.08 


405.5 


14.77 


5 X 10-2 


325.4 


345.2 


6.08 


330.0 


1.41 


5 X 10-^ 


252.9 


269.8 


6.68 


254.6 


0.67 




Bottom 


Oc 


(7exact (nb) 


a^w only log (nb) 


(%) 


CTww (nb) 


(%) 


5 X 10-1 


3.110 


3.404 


9.45 


3.277 


5.37 


5 X 10-2 


2.599 


2.732 


5.12 


2.605 


0.23 


5 X IC)-'"' 


1.931 


2.060 


6.(38 


1.933 


0.10 



Table 1: Total cross sections for charm and bottom production at HERA {E = 30 
GeV, = 314 GeV) for various angular cuts. The percentual differences between 
the exact results and the two approximate ones are also shown. Statistical errors are 
smaller that the last digit. 





Charm 


Oc 


(7exact (nb) 


CTww only log (nb) 


(%) 


CTww (nb) 


(%) 


5 X 10-1 


660.8 


866.4 


31.11 


839.0 


26.97 


5 X 10-2 


650.7 


732.7 


12.60 


705.3 


8.39 


5 X 10-^ 


549.6 


599.0 


8.99 


571.6 


4.00 




Bottom 




O-exact (nb) 


a^w only log (nb) 


(%) 


o-ivw (nb) 


(%) 


5 X 10-1 


8.585 


10.20 


18.81 


9.866 


14.92 


5 X 10-2 


7.997 


8.492 


6.19 


8.154 


1.96 


5 X 10-3 


6.360 


6.780 


6.60 


6.441 


1.27 



Table 2: Total cross sections for charm and bottom production ai E = 100 GeV and 
a/S* = 1 TeV for various angular cuts. The percentual differences between the exact 
results and the two approximate ones are also shown. Statistical errors are smaller 
that the last digit. 
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tion (|3l| ) gives 

g^.^ = -2k-q + AM^ = -2k ■ py + AM^ ~ -(ys - AM^). (32) 



The integration limits become 



-(ys-AM'') <g2<_!!^. (33) 
V ^ 1 

In this case, however, we cannot simply take = in Wi{q'^,k ■ q), as we did 
in the previous case. In fact, processes involving the absorption of a virtual photon 
by a hadron target are characterized by a dimensional parameter A, defined in such 
a way that the virtual photon cross section is nearly equal to the photoproduction 
cross section for < A^, while it falls rapidly to zero for virtualities above this 
value. In particular, for the production of heavy objects, one can assume A = ^M, 
where ^ is a dimensionless parameter of order 1, since the characteristic scale of the 
process is given by the mass of the heavy object. For this reason, with quite a drastic 
approximation, instead of eq. (^1]) we may assume 

W,{q^k-q) 



Ak ■ q 



a,,{q,k)e[q' + eM'), (34) 



which leads to 



Op 



' Ztt 
where the quantity 



2"^eZ/ I 7^ 1 + ' ' log 



Qww rii?y ) y mly'^ 



(35) 



Ql^ = mm{ys-AM\eM^) (36) 

is usually referred to as the effective Weizsacker- Williams scale. 

In practice, ys — AM'^ will also be of order M^, because of the competing effects of 
the large- a; suppression of the gluon density of the proton (which favours the threshold 
region) and dynamical threshold suppression. We therefore expect that in general the 
approximation 

Ql^ = (37) 

will be appropriate. We see that, at least for the production of heavy objects, the ef- 
fective Weizsacker- Williams scale is determined by dynamical rather than kinematical 
considerations. In the presence of an angular cut, dynamical arguments are ineffec- 
tive, because for small angles (which are of interest in photoproduction processes) the 
minimum value of the virtuality, eq. (p6D, is well below ^^M^ in absolute value. 
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The parameter ^ is essentially unconstrained (except that we expect it to be 
of order 1). One can therefore absorb the non-logarithmic term in eq. ( pSj ) by an 
appropriate redefinition of ^. Consider for example eq. ( p5D with only the most 
singular terms in y taken into account. We can get rid of the non-logarithmic term 
by replacing ^ with where ^' is determined by the condition 

2 2, i^M^ 2, i'^M^ , , 

- - + -log^-^ = -log^-^, (38) 

or 

^'-0.6^. (39) 

The numerical results for total cross sections at HERA, both for charm and bottom 
production, are presented in table ^ In the first row we show the result of the 
exact leading-order calculation. We see that the use of the usual Weizsacker- Williams 
approximation of eq. (P, shown in the second row, leads to a completely wrong result. 
The following two entries are obtained using eq. ( P5| ) with Q\,^ = ys—AM^. The effect 
of the inclusion of the non-logarithmic term is apparent, but still the approximation is 
not satisfactory. The remaining entries are obtained by implementing the dynamical 
considerations made above. We performed the calculation using only the logarithmic 
term in the Weizsacker- Williams distribution, and choosing different values for ^, 
namely ^ = 1/2,1,2. The choice ^ = 1 gives a very good agreement between the 
exact and the approximated result. We then tested eq. ( pPj ) by including the non- 
logarithmic term in the Weizsacker- Williams function and taking Qww = M/0.6. 
Finally, the comparison between the entries in the fifth and the last rows is a test of 
eq. (13). 
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Charm 


Bottom 




cr (nb) 

V / 


(%) 

V / 


a (nb) 


(%) 


Fjxart rpsiilt 


353.4 




3.135 




- 4M^ only log 


478.4 


35.37 


4.090 


30.46 


— 4M^ only log 


383.1 


8.40 


3.392 


8.20 


Vys - 4M^ 


367.9 


4.10 


3.265 


4.15 


M only log 


349.9 


0.99 


3.137 


0.06 


M/2 only log 


327.1 


7.44 


2.934 


6.41 


2M only log 


372.6 


5.43 


3.339 


6.51 


M/0.6 


351.4 


0.57 


3.158 


0.73 


min{\/ys — 4M^, M) only log 


349.0 


1.24 


3.128 


0.22 



Table 3: Total cross sections for heavy-flavour production at HERA, for different 
choices of the Weizsacker- Williams distribution. The percentual differences between 
the exact results and the approximate ones are also shown. 
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